On the scenario for the glass transition 
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In this letter we study a lattice glass system that undergoes a glass transition. When we approach 
the glass transition we find both a divergence of a point to set correlation length and a vanishing of 
the thermodynamic potential. These findings are similar to the predictions coming from mean-field 
(replica) theory, but they differ from these predictions in some details: they underline the need of 
a better theoretical understanding of the glass transition. 
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Second order phase transitions are usually character- 
ized by the existence of an equilibrium correlation length 
that is divergent at the phase transition point. It is usu- 
ally stated that such a correlation length is absent for the 
glass transition. However it was noticed long time that 
in the framework of the Adams Gibbs theory [TJ [H [3] 
the slow dynamics is dominated by large regions that 
move in a cooperative way. This observation has been 
formalized by introducing the appropriate correlations 
functions that can be used to define the dynamical cor- 
relation length [H |51 |SJ [HI [Hj ■ A precise theoretical anal- 
ysis has been done both in the framework the replica 
approach [5J [7J [TOJ EJ EH] and in the framework of the 
mode-coupling theory 13J . 

It was later realized that we can define a static correla- 
tion length that is divergent at the transition temperature 
[Til EES EES EES EES- It is related to the point to set cor- 
relations: its definition is the following. We start from a 
large system at equilibrium: we froze it and we allow only 
a region of size R to move. We ask if the equilibrium con- 
figuration satisfying this constraint is very far from the 
original configuration. It is natural to suppose that if 
R << £d> where £d is the typical size of the dynamically 
rearranging regions, the system is blocked. The opposite 
behaviour should be present for R » We can define 
an equilibrium correlation length £ s as the minimal size 
of a non blocked region: the previous arguments suggest 
that we should have £ s oc £<j . It has been shown [11] that 
in the mean field approach £ s diverges near the transition 
point. While real experiments to measure £ s are not yet 
feasible, numerical experiments have been done showing 
evidence of a large increase of the £ s near the thermody- 
namic transition point [15j. 

In this letter I will study a different correlation length 
£ w . The definition is similar to £ s , but it is technically 
different |16j : we froze the system on a wall (a two dimen- 
sional slice) and we measure how the final configuration 
is different from the initial one (the geometry is similar to 
the one used in |17j for the dynamics) . More precisely we 
introduce an overlap q that indicates the similarity of two 
configurations: q = 1 for identical configuration while q 
takes a small value for unrelated configuration. In this 
way, if we froze the region where < x < 1, x being one 




2 4 6 8 10 



FIG. 1: The data for q(x) - q bM at fj, = 8.5 (p = 0.8097) as 
function of I on a cubic lattice of side L = 22. The value of 
the correlation length is £ TO = 1.32 ± 0.03. 

of the coordinates, we can define the average value of q at 
given x, i.e. q(x). The wall correlation length is defined 

by 

q(x) sa aexp(-x/£ w ) + q bu ik (1) 

at large x. The correlation length £ w is logically different 
from £ s and its behaviour in mean field theory is not clear. 

In this letter we present evidence that the correlation 
length £ w is divergent when we approach the phase tran- 
sition in a lattice glass model (a similar conclusion has 
been reached for an off lattice fluid [H]). We will relate 
the divergence of £ w to the development of flat direction 
for an appropriate thermodynamic potential [18 . 

The model we consider is a simple lattice gas model 
(the CTCC model [H]), that does not easily crystallize in 
simulations (the cell crystal contains 7 3 = 343 particles 
and it is quite difficult to form). Each site of a three 
dimensional cubic lattice may be empty or occupied. If 
the site is occupied we have to specify the orientation 
of the particle. The particles may point in each of the 
six directions that correspond to elementary moves on 
the lattice. The system satisfies the constraints: 1) each 
lattice point may occupied by only one kind of particles, 
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2) if the site i is occupied by a particle that points in the 
direction d, the site i + d is empty. 

We introduce a chemical potential fi. In the limit 
/i — > oo we reach the maximum density state, i.e. a 
crystal with periodicity 7 and density 6/7=0.8571. In 
simulations on large lattices, starting from the disorder 
phase, there is no sign of crystallization up to the longest 
runs. One observes in the dynamics a very nice slowing 
down, with all the standard stigmata of a glass transition, 
i.e. mode-coupling behaviour of the correlation function 
(with both the a and the f3 region well exposed) , a strong 
increase of the dynamical susceptibilities and a simulta- 
neous decrease of the the diffusion constant D pjjj. This 
data on the diffusion constant do not have a particular 
nice behaviour as function of the chemical potential fi. 
They are much more regular as function of the density. 

We will show that there are statically defined quan- 
tities that display a critical behaviour: they should be 
counterparts of these dynamical quantities. Let as start 
with a precise definition of £ w . We consider a large lat- 
tice of size L: the value of L will not enter in the analysis 
but consistency implies that L » £ w . At given chemical 
potential /i we take a thermalized configuration a. We 
now consider a system where at x = the configurations 
t satisfy the constraint 

T{0,y,z) = a(Q,y,z) . (2) 

In other words the particles on the wall (i.e. at x = 0) 
cannot move and the particles outside the wall cannot 
penetrate inside. 
We define 

q„{x,y,z) = (q a ,T(x,y,z)) T , q(x,y,z) = (q a (x, y, z)) a .(3) 

Here (•)„• and (•),- denote respectively the statistical aver- 
age over the variables a and r. We perform a double sta- 
tistical average, first on the r and later on the a. Trans- 
lational invariance implies that q(x, y, z) does not depend 
on y and z and it will be denoted by q(x). This construc- 
tion is similar to the one used in [TTJ Q21 H3J HH H] for 
defining £ s . Here the constraint eq. Q is imposed on a 
slice an most of the system is free; in the definition of £ s 
the same constraint is imposed on the whole system and 
only a small region is left free. 

The function q{x) tell us how much two equilibrium 
configurations are different if they coincide at distance x 
in one direction. Which should be the behaviour of q(x) 
near the thermodynamic glass transition? In the mosaic 
picture the function q{x) should be strongly different from 
its asymptotic value at large x in a region that becomes 
larger and larger when we approach the glass transition. 
The shape of the function q(x) near the wall is an inter- 
esting problem that we do not address here. We focalize 
our interest on the behaviour at large x. We find that for 
x > 2 the function q{x) can be remarkably well fitted by 
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FIG. 2: The correlation length £ w in the region of p ranging 
from 1 to 9.25 as function of p — p c , with p c = 0.826. 

an exponential, i.e. 

q{x) - qbuik = aexp(-x/£ w ) + c (4) 

where c is a small constant, i.e. O(10 -3 ) (the previous 
hand-waving argument suggest that £ w should diverges 
at the thermodynamic phase transition were £ s diverges). 
An example of the fit is shown in fig.Q, where the data 
are taken at fi = 8.5 (p = 0.810) on a cubic lattice of side 
L = 22. The value of the correlation length is 1.32 ±0.03. 

We have measured this correlation (or penetration) 
length £ w in the region of /1 ranging from 1 to 9.25 on 
lattices of various side (up to L — 28). The correlation 
£ w is very well fitted by a power low behaviour 

Up) = a(pc-pV , (5) 

with p c = 0.826±0.002 and v = 0.48±0.02 (see fig. Q). 
The value of the exponent v is very near to the simple 
value 1/2. 

In the study of phase transitions one usually introduces 
an order parameter and an associated thermodynamic po- 
tential. If the order parameter is a continuos function of 
the temperature, at the phase transition point the poten- 
tial is zero only in one point, while at lower temperatures 
it is different from zero in a region whose side goes to 
zero at the critical point. For a transition with a discon- 
tinuous order parameter, at the phase transition point 
the potential is zero when the order parameter stays in 
a given interval and a similar behaviour is present also 
below the phase transition point. In the infinite volume 
limit the potential must be convex function: when the 
potential is not convex in a mean field approximation, 
Maxwell construction enforces the convexity. 

The glass transition may be characterized by the exis- 
tence of regions of configuration space where (below the 
transition) the system remains trapped for an infinitely 
large time. The overlap q is the putative order parameter 
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FIG. 3: The same data as fig |2|: here we plot £, w {p) 2 versus 
p in a linear scale. 
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FIG. 4: The derivative of the potential W' p (q) at various den- 
sities (0. 749, 0.771, 0.789, 0.803, 0.810) on large systems of 
side L — 40. 



of the glass transition. One can sharpen the physical pic- 
ture by introducing a potential W(q) defined as follows 
[IH1 [20J • We consider an equilibrium configuration a. 
We call P a (q) the probability that an other configuration 
r has an overlap q a T = q. We define 



W(q) 



lim 



v 



(6) 



With probability one when the volume V goes to infin- 
ity, the potential W(q) does not depend on the reference 
configuration a. In other words P a (q) ~ exp(— VW(q)). 
By construction W(qb u ik) = and the vanishing of the 
potential W(q) for more than one q-value is the distin- 
tive characteristic of replica symmetry breaking (it should 
happen below the thermodynamic glass transition). 

There are many possible techniques that can be used 
to compute the potential W(q). One of the best is the 
introduction of a tethering potential, as described in ffi4"] . 



Here we use a simpler method |25j : we constrain the r 
configuration to stay in the region q a T > q*. With this 
constrain the resulting probability distribution of q aiT is 
proportional to 9{q—q*) exp(— VW(q— q*)). In the region 
where q* > qbuik, we have that 



V(q^ T - q*) = 



V dq 



= l/W'{q*) 



(7) 



Using a wild metaphor we could say that W'(q) corre- 
spond to the force that pushes the system away from the 
region where q is high. It is interesting to note that in 
the Maxwell construction W'(q) is constant in the region 
where the naive mean field prediction violates convexity. 

In the replica approach to glasses the potential W(q) 
plays a central role [B]: in the mean field approach at 
high temperature the potential is convex. By decreasing 
the temperature firstly the convexity is lost; at a lower 
temperature (i.e. at the dynamic or mode-coupling tran- 
sition) a secondary local minimum (q m ) appears. The 
height of the minimum (i.e W(q m )) is the configurational 
complexity S. By decreasing the temperature the value 
of W(q) at the secondary minimum decrease up to the 
point (the Kauzmann transition) where W(q m ) = S = 0. 
Below this temperature we are in the glass phase: two 
minima (with W{q) — 0) are present and the replica sym- 
metry is broken. This scenario is valid in infinite range 
model where very detailed computation can be done. 

It is clear that the previous picture must be strongly 
modified in finite dimensions. We have already remarked 
the any thermodynamic potential in the infinite volume 
limit must be convex. A consistent scenario could be pro- 
vided by the Maxwell construction, unfortunately there 
are not many numerical result on the behaviour of the 
potential in finite dimensional models [6j. 

In this letter we compute the the potential W p (q) in the 
CTCC lattice gas model and we find that it has some sur- 
prising properties. The derivative of the potential W'{q) 
at various densities (from 0.749 to 0.810) is shown in figs. 
Q and ([5]). A few remarks are in order: 

• We cannot identify regions where W'Jq) decreases 
by increasing q: convexity is always strictly satis- 
fied. 
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are no flat regions {W'Jq) — const) that 
could be the remnant of the Maxwell construction. 
Therefore no dynamical (mode coupling) transition 
can be identified. 

• The potentials at different densities have similar 
shape. 

• The potential strongly decrease in nearly the whole 
q interval when the density increases. 

The last two statements can be done more sharply if 
we rescale the potential with the square of the correlation 
length 



n'(q) = Up) 2 W'M ■ 



(8) 
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FIG. 5: The same data of fig. Q for derivative of the potential 
W' p (q) multiplied by £™(p) 2 , i.e. Q' p (q). 

The data for Cl' p (q) are shown in fig. ([51 and they are 
weakly dependent on the density. If £ w (p) goes to infinity 
at a critical value of p the potential W(q) goes to zero in 
the whole interval from q m up to a value of q of near to 
0.9. 

The results are surprising (at least to me) for the fol- 
lowing reasons: 

• The potential W(q) cannot be used to identify the 
(dynamic) mode coupling transition: the configura- 
tional complexity cannot be defined. However it is 
possible that in this model the dynamic transition 
is very near to the static one, so that the two transi- 
tions cannot be easily disentangled: is this case one 
may observing overlapping features. One should 
also remember that W(l) is the entropy. If at the 
Kauzmann transition the entropy takes a very small 
value, the potential W(q) should be nearly vanish- 
ing in large interval of q. 

• The extrapolation of these results implies that a 
static Kauzmann transition is present and the be- 
haviour is rather simple. The exponent v is very 
near to the simple value 1/2 and the thermody- 
namic Wp(q) potential scales naively without any 
need of introducing extra power corrections. 

Similar investigations are presently done also for off- 
lattice fluids [T5J [55] and it would be very interesting to 
compare the results. A very important question is the 
theoretical derivation of these results. 
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